Quantum Spin Hall Effect 
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The quantum Hall liquid is a novel state of matter with profound emergent properties such as 
fractional charge and statistics. Existence of the quantum Hall effect requires breaking of the time 
reversal symmetry caused by an external magnetic field. In this work, we predict a quantized 
spin Hall effect in the absence of any magnetic field, where the intrinsic spin Hall conductance is 
quantized in units of 2^-. The degenerate quantum Landau levels are created by the spin-orbit 
coupling in conventional semiconductors in the presence of a strain gradient. This new state of 
matter has many profound correlated properties described by a topological field theory. 



Recently, the intrinsic spin Hall effect has been the- 
oretically predicted for semiconductors with spin-orbit 
coupled band structures Q, Q- The spin Hall current is 
induced by the external electric field according to the 
equation 
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where jj is the spin current of the i-th component of the 
spin along the direction j, Ek is the electric field and e^fc 
is the totally antisymmetric tensor in three dimensions. 
The spin Hall effect has recently been detected in two 
different experiments^^], and there is strong indication 
that at least one of them is in the intrinsic regime 0. 
Because both the electric field and the spin current are 
even under time reversal, the spin current could be dis- 
sipationless, and the value of o~ s could be independent 
of the scattering rates. This is in sharp contrast with 
the extrinsic spin Hall effect, where the effect arises only 
from the Mott scattering from the impurity atoms 0. 

The independence of the intrinsic spin Hall conduc- 
tance a s on the impurity scattering rate naturally raises 
the question whether it can be quantized under certain 
conditions, similar to the quantized charge Hall effect. 
We start off our analysis with a question: Can we have 
Landau Level (LL) -like behavior in the absence of a mag- 
netic field 0? The quantum Landau levels arise physi- 
cally from a velocity dependent force, namely the Lorentz 
force, which contributes a term proportional to A ■ p in 
the Hamiltonian. Here p is the particle momentum and 
A is the vector potential, which in the symmetric gauge 
is given by A = — x,0). In this case, the velocity 
dependent term in the Hamiltonian is proportional to 
B(xp y -yp x ). 

In condensed matter systems, the only other ubiqui- 
tous velocity dependent force besides the Lorentz force 
is the spin-orbit coupling force, which contributes a term 
proportional to (px E) ■ a in the Hamiltonian. Here E is 
the electric field, and a is the Pauli spin matrix. Unlike 
the magnetic field, the presence of an electric field does 
not break the time reversal symmetry. If we consider the 
particle momentum confined in a two dimensional geom- 
etry, say the xy plane, and the electric field direction con- 
fined in the xy plane as well, only the z component of the 
spin enters the Hamiltonian. Furthermore, if the electric 



field E is not constant but is proportional to the radial 
coordinate r, as it would be, for example, in the inte- 
rior of a uniformly charged cylinder E ~ E(x, y, 0), then 
the spin-orbit coupling term in the Hamiltonian takes the 
form Eo- z (xp y — yp x )- We see that this system behaves 
in such a way as if particles with opposite spins experi- 
ence the opposite "effective" orbital magnetic fields, and 
a Landau level structure should appear for each spin ori- 
entations. 

However, such an electric field configuration is not easy 
to realize. Fortunately, the scenario previously described 
is realizable in zinc-blende semiconductors such as GaAs, 
where the shear strain gradients can play a similar role. 
Zinc-blende semiconductors have the point-group sym- 
metry Td which is half of the cubic-symmetry group Oh, 
and does not contain inversion as one of its symmetries. 
Under the Td point group, the cubic harmonics xyz trans- 
form like the identity, and off-diagonal symmetric tensors 
(xy + yx, etc.) transform in the same way as vectors on 
the other direction (z, etc), and represent basis func- 
tions for the Ti representation of the group. Specifically, 



strain is a symmetric tensor 



and its off-diagonal 



(shear) components are, for the purpose of writing down 
a spin-orbit coupling Hamiltonian, equivalent to an elec- 
tric field in the remaining direction: 
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The Hamiltonian for the conduction band of bulk zinc- 
blende semiconductors under strain is hence the anal- 
ogous to the spin-orbit coupling term (v X E) ■ a. In 
addition, we have the usual kinetic p 2 term and a trace 
of the strain term tre = e xx + e yy + e zz , both of which 
transform as the identity under T^: 
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For GaAs, the constant ^ = 8x 10 5 m/s This Hamil- 
tonian is not new and was previously written down in Ref. 
El EH E3 but the analogy with the electric field and 
its derivation from a Lorentz force is suggestive enough to 
warrant repetition. There is yet another term allowed by 
group theory in the Hamiltonian but this is higher 



2 



order in perturbation theory and hence not of primary 
importance. 

Let us now presume a strain configuration in which 
£xy = but e xz has a constant gradient along the y di- 
rection while e yz has a constant gradient along the x di- 
rection. This case then mimics the situation of the elec- 
tric field inside a uniformly charged cylinder discussed 
above, as e xz (<-» E y ) = gy and e yz (^ E x ) = gx, g being 
the magnitude of the strain gradient. With this strain 
configuration and in a symmetric quantum well in the 
xy plane, which we approximate as being parabolic, the 
above Hamiltonian becomes: 



H = 
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We first solve this Hamiltonian and come back to the 
experimental realization of the strain architecture in 
the later stages of the paper. We make the coordi- 



nate change x — > (2mD) 



-1/4 



x, y 



{2mD)- 1 / i y and 



R = 5xV^ 5 ' R = 2 ° r D = D ° = 
a special point, where the Hamiltonian can be writ- 
ten as a complete square, namely H — ^(p — eA<j z ) 2 
with A — ™£f? 9 (y, —x, 0). At this point, our Hamilto- 
nian is exactly equivalent to the usual Hamiltonian of 
a charged particle in an uniform magnetic field, where 
the two different spin directions experience the oppo- 
site directions of the "effective" magnetic field. Any 
generic confining potential V(x, y) can be written as 
Dq{x 2 + y 2 ) + AV(x,y), where the first term completes 
the square for the Hamiltonian, and the second term 
AV(x,y) = V(x,y) — D (x 2 + y 2 ) describes the addi- 
tional static potential within the Landau levels. Since 
[H, <t z ] = we can use the spin on the z direction to char- 
acterize the states. In the new coordinates, the Hamilto- 
nian takes the form: 
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The -ffj. i is the Hamiltonian for the down and up spin er 2 
respectively. Working in complex-coordinate formalism 



and choosing z = x - 
lowering operators: 

a = d z 
b = d z 



■ iy we obtain two sets of raising and 
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in terms of which the Hamiltonian decouples: 



(l T |)oot + (l±|)66t + l 




(6) 



(7) 



The eigenstates of this system are harmonic oscil- 
lators \m,n) = (at) m (6t)«|o,0) of energy E$ n = 



^[(lT#)m+(l±#)n+l] 
the case of R = 2 where there is no additional static 
potential within the Landau level. 

For the spin up electron, the vicinity of R w 2 is charac- 
terized by the Hamiltonian — ^^g{2aa) + 1) with the 
LLL wave function ^(z) — ^= exp( ~y ). a is the op- 
erator moving between different Landau levels, while b is 
the operator moving between different degenerate angu- 



We shall focus on 



lar momentum states within the same LL: L, 



btf- 



-aa> 



Lzty^iz) = n<jy n (z). The wave function, besides the con- 
fining factor, is holomorphic in z, as expected. These up 
spin electrons are the chiral, and their charge conduc- 
tance is quantized in units of e 2 /h. 

For the spin down electron, the situation is exactly 
the opposite. The vicinity of R s» 2 is characterized 
by the Hamiltonian H L = \^g(2btf + 1) with the LLL 

wave function 4>ln{z) = % ^ , exp( ~y ). b is the oper- 
ator moving between different Landau levels, while a is 
the operator between different degenerate angular mo- 
mentum states within the same LL: L z — btf — aa^ , 
L z 4>l n (z) = —m4>l n (z). The wave function, besides the 
confining factor, is anti-holomorphic in z. These down 
spin electrons are anti-chiral, and their charge conduc- 
tance is quantized in units of —e 2 /h. 



u z ~xy 




FIG. 1: Spin up and down electrons have opposite chiral- 
ity as they feel the opposite spin-orbit coupling force. Total 
charge conductance vanishes but spin conductance is quan- 
tized. The inset shows the lattice displacement leading to the 
strain configuration. 

The picture that now emerges is the following: our 
system is equivalent to a bilayer system; in one of the 
layers we have spin down electrons in the presence of a 
down-magnetic field whereas in the other layer we have 
spin up electrons in the presence of an up-magnetic field. 
These two layers are placed together. The spin up elec- 
trons have positive charge Hall conductance while the 



3 



spin down electrons have negative charge Hall conduc- 
tance. As such, the charge Hall conductance of the whole 
system vanishes. The time reversal symmetry reverses 
the directions of the "effective" orbital magnetic fields, 
but interchanges the layers at the same time. However, 
the spin Hall conductance remains finite, as the chiral 
states are spin-up while the anti-chiral states are spin- 
down, as shown in FigurepQ. The spin Hall conductance 
is hence quantized in units of 2^-^ = 2j-. Since an 
electron with charge e also carries spin H/2, a factor of 
^ is used to convert charge conductance into the spin 
conductance. 

While it is hard to experimentally measure the spin 
Hall effect, and supposedly even harder to measure the 
quantum version of it, the measurement of the charge 
quantum Hall effect has become relatively common. In 
our system, however, the charge Hall conductance a xy 
vanishes by symmetry. However, we can use our physi- 
cal analogy of the two layer system placed together. In 
each of the layers we have a charge Quantum Hall effect 
at the same time (since the filling is equal in both lay- 
ers), but with opposite Hall conductance. However, when 
on the plateau, the longitudinal conductance a xx also 
vanishes (a xx — 0) separately for the spin-up and spin- 
down electrons, and hence vanishes for the whole system. 
Of course, between plateaus it will have non-zero spikes 
(narrow regions). This is the easiest detectable feature 
of the new state, as the measurement is entirely electric. 
Other experiments on the new state could involve the in- 
jection of spin-polarized edge states, which would acquire 
different chirality depending on the initial spin direction. 

We now discuss the realization of a strain gradient of 
the specific form proposed in this paper. The strain 
tensor is related to the displacement of lattice atoms 
from their equilibrium position Ui in the familiar way 
€ij = (dui/dxj + duj/dxi)/2. Our strain configuration 
is e xx = £ yy = t Z z = f-xy = as well as the strain gra- 
dients e zx — gy and e yx — gx. Having the diagonal 
strain components non-zero will not change the physics 
as they add only a chemical potential term to the Hamil- 
tonian. The above strain configuration corresponds to a 
displacement of atoms from their equilibrium positions 
of the form u — (0,0,2gxy). This can be possibly re- 
alized by pulverizing GaAs on a substrate in MBE at a 
rate which is a function of the position of the pulveriz- 
ing beam on the substrate. The GaAs pulverization rate 
should vary as xy ~ r 2 sin(2</>), where r is the distance 
from one of the corners of the sample where the GaAs 
depositing was started. Conversely, we can keep the pul- 
verizing beam fixed at some r and rotate the sample with 
an angle-dependent angular velocity of the form sin(20). 
We then move to the next incremental distance r, in- 
crease the beam rate as r 2 and again start rotating the 
substrate as the depositing procedure is underway. 

The strain architecture we have proposed to realize the 
Quantum Spin Hall effect is by no means unique. In the 



present case, we have re-created the so-called symmet- 
ric gauge in magnetic-field language, but, with different 
strain architectures, one can create the Landau gauge 
Hamiltonian and indeed many other gauges. The Landau 
gauge Hamiltonian is maybe the easiest to realize in an 
experimental situation, by growing the quantum well in 
the [110] direction. This situation creates an off-diagonal 
strain e xy — \SaaT, and e xz = e yz — where T is the 
lattice mismatch (or impurity concentration), S44 is a ma- 
terial constant and x, y, z are the cubic axes. The spin- 
orbit part of the Hamiltonian is now ^-e xy (p x a y —p y cr x ). 
However, since the growth direction of the well is [110] 
we must make a coordinate transformation to the x',y', z' 
coordinates of the quantum well {x' , y' are the new coor- 
dinates in the plane of the well, whereas z' is the growth 
coordinate, perpendicular to the well and identical to the 
[110] direction in cubic axes). The coordinate transfor- 
mation reads: x' = -j=(x — y) , y' = —z, z' — -^(x + y), 
and the momentum along z' is quantized. We now vary 
the impurity concentration T (or vary the speed at which 
we deposit the layers) linearly on the y' direction of the 
quantum well so that e xy = gy' where g is strain gradient, 
linearly proportional to the gradient in T. In the new co- 
ordinates and for this strain geometry, the Hamiltonian 
reads: 

H=^ + ^gy' Px >cr z >+Dy 2 (8) 
2m a 

where we have added a confining potential. At the suit- 
able match between D and <?, this is the Landau-gauge 
Hamiltonian. One can also replace the soft-wall condition 
(the Harmonic potential) by hard-wall boundary condi- 
tion. 

We now estimate the Landau Level gap and the strain 
gradient needed for such an effect, as well as the strength 
of the confining potential. In the case R w 2 the en- 
ergy difference between Landau levels is AE L andau — 
2 x ^\^h9 = ^39- For a gap of lmK : we hence need 
a strain gradient or 1% over 60/ito. Such a strain gra- 
dient is easily realizable experimentally, but one would 
probably want to increase the gap to lOmK or more, 
for which a strain gradient of 1% over 6/im or larger is 
desirable. Such strain gradients have been realized ex- 
perimentally, however, not exactly in the configuration 
proposed here [Til Il5| . The strength of the confining po- 
tential is in this case D — 10~ 15 7V/m which corresponds 
roughly to an electric field of 1 V/m for a sample of 60/xm. 
In systems with higher spin-orbit coupling, C3 would be 
larger, and the strain gradient field would create a larger 
gap between the Landau levels. 

We now turn to the question of the many-body wave 
function in the presence of interactions. For our system 
this is very suggestive, as the wave function incorporates 
both holomorphic and anti-holomorphic coordinates, by 
contrast to the pure holomorphic Laughlin states. Let 
the up-spin coordinates be Zi while the down-spin coordi- 
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nates be the Wi . Zi enter in holomorphic form in the wave 
function whereas Wi enter anti-holomorphically. While if 
the spin-up and spin-down electrons would lie in sepa- 
rate bi-layers the many-body wave function would be just 
IW* - Zj) m Tlk<M ~ «;?re-4C 4 .-?+E*^J), 
where m is an odd integer. Since the particles in our 
state reside in the same quantum well and may possibly 
experience the additional interaction between the differ- 
ent spin states, a more appropriate wave function is: 

U rs (z r - l( ,J)"e4E. 2 ' 2 .'+^ ,t ^) (9) 

The above wave function is symmetric upon the inter- 
change z <-> w* reflecting the spin-| chiral - spin- J. anti- 
chiral symmetry. This wave function is of course anal- 
ogous to the Halperin's wave function of two different 
spin states . The key difference is that the two differ- 
ent spin states here experience the opposite directions of 
magnetic fields. 

Many profound topological properties of the quantum 
Hall effect are captured by the Chern-Simons-Landau- 
Ginzburg theory [13. While the usual spin orbit coupling 
for spin- 1/2 systems is T- invariant but P-breaking, our 
spin-orbit coupling is also P-invariant due to the strain 
gradient. The low energy field theory of the spin Hall 
liquid is hence a double Chern-Simons theory with the 
action: 

S = £ / ^"apdvap - £ J e^c^cp (10) 

where the and c M fields are associated with the left 
and right movers of our theory while v is the filling fac- 
tor. The fractional charge and statistics of the quasi- 
particle follow easily from this Chern-Simons action. Es- 
sentially, the two Chern-Simons terms have the same fill- 
ing factor v and hence the Hilbert space is not the tensor 
product of any two algebras, but of two identical ones. 
This is a mathematical statement of the fact that one 
can insert an up or down electron in the system with the 
same probability. Such special theories avoids the chiral 
anomaly |18| and their Berry phases have been recently 
proposed as preliminary examples of topological quan- 
tum computation [l^. It is refreshing to see that such 
abstract mathematical models can be realized in conven- 
tional semiconductors. 

A similar situation of Landau levels without magnetic 
field arises in rotating BECs where the mean field Hamil- 
tonian is similar to either or Hi . In the limit of rapid 
rotation, the condensate expands and becomes effectively 
two-dimensional. The L z term is induced by the rotation 
vector [IS!. The LLL behavior is achieved when the 
rotation frequency reaches a specific value analogous to 
the case R « 2 in our Hamiltonian. In the BEC litera- 
ture this is the so called mean-field quantum Hall limit 



and the ground state wave function is Laughlin-type. In 
contrast to our case, the theory is still T breaking, the 
magnetic field is replaced by a rotation axial vector, and 
the lowest Landau level is chiral. 

In conclusion we predict a new state of matter where a 
quantum spin Hall liquid is formed in conventional semi- 
conductors with spin-orbit coupling. The quantum Lan- 
dau levels are caused by the gradient of strain field, rather 
than the magnetic field. The new quantum spin Hall liq- 
uid state shares many emergent properties similar to the 
charge quantum Hall effect, however, unlike the charge 
quantum Hall effect, our system does not violate time 
reversal symmetry. 
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